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Abstract. We present an XMM observation of the distant galaxy cluster CL 0939-1-4713. The X-ray image 
shows pronounced substructure, with two main subclusters which have even some internal structure. This is 
an indication that the cluster is a dynamically young system. This conclusion is supported by the temperature 
distribution: a hot region is found between the two main subclusters indicating that they are at the beginning 
of a major merger, and that they will probably collide in a few hundreds of Myr. The intra-cluster gas of 
CL 0939-1-4713 shows inhomogeneities in the metal distribution, with the optically richer subcluster having a 
higher metallicity. 
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1. Introduction 

The galaxy cluster CL 0939-^4713 (Abell 851) at a redshift 
z = 0.41 is an exceptionally rich cluster. It was extensively 
observed in morphological and spectroscopic studies of 
the galaxy population. These stud ies are based on ground- 
based o ptical and near-infrared dDressler fc Gunnlll992l: 
Belloni eTal] Il995t iFuku gita et alJ Il995t IStanford et all 
19951: iBelloni fc Roseil Il996l: ~E)ressler et al.l Il999t 



Poggia nti et all Il999t live et alJ I2OOOI: iMartin et all 
_2000') and HST ob s ervations jPressler et al.lll99lll994t 
Smail et alJ Il999t IZiegler et alJ Il999t iFerrer as fc Silk 
as well as obs ervations in UV llBuson et al.ii200Q|) 
and sub-millimeter ijCowie et al.l 1200^ wavelengths. A 
high number of star-burst and post-starburst galaxies was 
found in this distant cluster. To study the evolution of 
galaxies this number was compared with the correspond- 
ing number in nearby clusters. In the surroundings of the 
cluster s a network of galaxy filaments and subclumps was 
found l|Kodama et aLlboOljl . 

The cluster was a lso the target of several gravitational 
lensing analyses llSeitz et alJ Il99(l iGeiger fc Schneide^ 
Il999t live etaJ] l200^^ aiming at the determination 
of the cluster mass distribution. The comparison be- 
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tween the masses det e rmine d by the weak lensing 
an alysis bv ISeitz et alJ lll996j) and the X-ray analysis 
by ISchindler et alJ l)l998() showed a discrepancy of a 
factor of 2-3. 

Previous X-ray observations by ROSAT and ASCA 
showed that the cluster has a very irregular structure 
indicatin g that the cluster is not relaxed but dynamically 
young ( S chindler fc WambsganssI Il996l: ISchindler et alJ 
|l998): this result is in good agreement wit h the c onclu- 
sion from larger scales by iKodama et al.l l|200l|) . The 
formation and evolution of clusters depends sensitively 
on cosmological para meters like the mean matter density 
in th e universe flm, ( Thomas et al.l Il998l: Ijenkins et alJ 
Il998t iBeisbart etaO boOlj) . Therefore it is important 
to determine the dynamical state of clusters at differ- 
ent redshifts, i.e. at different evolutionary states. The 
X-ray morphology alone is not the best indicator of the 
dynamical state, but it should be complemented with 
all other information available, e.g. the temperature 
map or the galaxy distribution. The combination of all 
findings gives a detailed picture of the state of a cluster. 
We do this analysis here for the moderately distant 
cluster CL 0939-1-4713, which can be used as a basis for 
comparisons at lower and higher redshifts. 
The capability of XMM to perform spatially resolved 
spectroscopy can be used also to determine the distri- 



2 



E. De Filippis et al.: XMM observation of CL 0939+4713 



bution of the metal abundances. Not only the overall 
value of metallicity but also its spatial distribution gives 
important indications on the metal enrichment processes 
and the cluster formation process. Early enrichment 
would point pre ferentially toward s galactic winds driven 
by supernovae (jPe Yound Il978|) . For this process one 
would expect a priori no metallicity variations, because 
there was enough time to distribute the metals uniformly 
throughout the cluster. However, numerical simula- 
tions of this mecha nism ga ve very d iscordant results. 
lAguirre et aD l(200l[) and .Metzler fc Evrard ( 1994, 1997) 
concluded that winds are very efficient in enriching the 
ICM, but the latter authors found a very steep radial 
metallicity gradient , which is not observed in any cluster. 
iMiirakami fc Bab^ill l(l999l) and ICnedinl l(l998li on the 
other hand found that galactic winds play only a minor 
role. 

At a later stage ram pressure stripping l|Gunn fc GottI 
Il972^ is expected to become increasingly important. If 
considerable amounts of metals are stripped recently 
a spatial variation - not necessarily radial gradients - 
should be observable. To test these effects it is particularly 
important to analyse clusters at different redshifts to see 
how the efficiency of the different enrichment mechanisms 
changes with time. 

Throughout this paper we use Hq = 50 km s~ Mpc^ 
and qo — 0.5; all errors are 90% confidence levels. 



2. Observation and data analysis 

2.1. Observation and data processing 

The XMM-Newton observation of the galaxy cluster 
CL 0939+4713 (ID 0106460101) was performed with the 
European Photon Imaging Camera (EPIC) on November 
6th, 2000 during the satellite's 167'^ revolution. It 
consists of one pointing obtained operating the EPIC 
cameras in the standard Full Frame Mode. For all the 
three EPIC cameras (MOSl, M0S2, and pn) the thin 
filter was used for a total exposure time of 51 ks. 
We used the SASv5.0.1 processing tasks emchain and 
epchain to generate calibrated event files from raw data. 
The XMM-Newton observations are strongly affected by 
high internal and external background levels. Efficient 
ways of subtracting these background compo nents 
were ex tensively discussed by IPratt et al.l l|2002l) and 
iLumb et al. (2002) . 
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Fig. 1. MOSl light curve in the energy range 10 — 12 keV. 
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Fig. 2. Effective exposure time versus the threshold in the 
count rate in the energy band 10 — 12 keV and 12 — 14 keV 
for MOSl & M0S2 and for pn, respectively. 



for pn) we excluded all the intervals of exposure time hav- 
ing a count rate higher than a certain threshold value. 



2.2. Flare rejection 

Low energy (< 1 MeV) protons from solar flares are 
stopped in the first microns of silicon of the EPIC CCDs. 
They give a signal identical to X-rays, which makes it 
impossible to recognize them via a pattern selection. We 
obtained a cleaner set of data by visually inspecting the 
light curves and applying an intensity filter. In the energy 
range 10-12 keV for MOSl & M0S2 cameras (12-14 keV 



The three light curves, one for each instrument, were 
analyzed separately in order to find the best threshold val- 
ues (see Fig.0. Different threshold values lead to different 
cuts in the exposure time (see Fig. ^ ; the selected values 
were therefore a compromise between maximizing the ex- 
posure time and minimizing the presence of events caused 
by fiares. 

The chosen threshold values are 15 counts -100 s^^ for 
MOSl and M0S2 and 25 counts -100 s^^ for pn; the 
threshold for the pn camera is higher because of the larger 
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Fig. 3. Left: Projection of the cluster surface brightness extracted from the rectangular region shown in the right 
panel. Right: X-ray image of the cluster. Superposed are the linearly spaced X-ray contours, the identification numbers 
of the point sources detected in the field (see § IH.4|I and the location of the two subclusters forming the cluster core: 
Ml and M2, correspondent to the two maxima in the X-ray emission. 



effective area of the instrument compared to MOS. The re- 
sulting effective exposure times for the three cameras are 

MOSl: 45800 s 
M0S2: 46300 s 
pn: 41800 s 

2.3. Vignetting correction 

To co rrect for vignetting t he method described in de- 
tail bv lArnaud et al.l l)2001a(l has been used. A weighting, 
given by the ratio of the effective area at the centre of 
the detector -^o,o(£'i)- to the effective area at the event 
position -Ax-^y^ (Ei)-, is assigned to each photon extracted 
in each spectrum. The vignetting corrected photon counts 
-C{Ei)- in each energy channel Ej {Ei — AEi/2 < Ej < 
E, + AEi/2) are 

j 

with the weights Wj being, as described above: 

Ao^ojEj) 
^ Axj,vj{Ei) 

which therefore all have values greater than unity. 
The errors in the spectra are recomputed taking into ac- 
count the vignetting correction; the variances in each en- 



ergy channel are therefore given by 

i 

The above method is easily generalized for imaging, as 
described bv FMaicrowicz et al.. (,2002,1 . 

2.4. Background subtraction 

In order to subtract both the cosmic an d the non cosmic 
X-ray background l|Arnaud et al.ll2002|) . the blank field 
observations, collected by D. Lumb were used. ^ These 
are homogeneous collections of high Galactic latitude 
observations where no bright source is present, for a total 
exposure time of about 400 ksec. These observations were 
all taken in a period of very low flare occurrence and 
almost exclusively with the use of a thin filter; bad pixels 
and sources were also removed. 

We applied to the background event files the same in- 
tensity filter applied to our observations; this yields final 
exposure times for the blank fields of MOSl: 373425 s, 
M0S2: 358904 s, and pn: 332646 s. Since the cosmic 
ray induced background might slightly change with time, 
the count rate in the 10.0 - 12.0 keV (12.0 - 14.0 keV) 

^ Event files can be downloaded via anony- 
mous ftp at xmm.vilspa.esa.es in the directory 
/ pub /ccf/constituents / extras /background. 
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energy band for MOSl & M0S2 (pn) cameras in the 
blank field was compared with the count rate in our 
observation (after flare rejection). The obtained normal- 
ization factors (0.995, 0.992 and 1.038 for MOSl, M0S2 
and pn, respectively) are then used to renormalize the 
blank field data. 

We then looked for a possible residual soft background 
due to the soft component (E < 1.5 keV) of the cosmic 
X-ray diffuse emission, which varies with the position in 
the sky. At the position in the sky of CL 0939-1-4713, 
no excess of soft cos mic background was observed in 
previous observations ijSnowden et al.lll997|) . 
Nevertheless, we performed a further check. We chose a 
region in our pointing where no emission from the cluster, 
or from other point sources, was observed. For all the 
three detectors, spectra were extracted from this region, 
together with background spectra from the blank fields 
by D. Lumb. No evident excess of flux, associated with 
the soft X-ray background, is observed in any of the three 
detectors. We then subtracted the blank field spectra 
from the ones extracted in our observation; the resulting 
spectra yield no significant signal and arc dominated by 
noise only. 

No additional correction for excess of soft background is 
therefore needed for CL 0939-^4713. 



3. Morphological analysis 

Figure El shows an image of CL 0939-1-4713 in the 
0.3 — 2.0 keV band. The image is an adaptively smoothed 
and exposure corrected mosaic of photon events from the 
three detectors (MOSl, M0S2 and pn). 
The cluster emission has an elliptical shape, centred at 
RA= 9^ 43™ 01.00% DEC= -h46° 59' 37.50" (J2000). 
The axis ratio of major to minor axis is 1.3. The position 
angle of the ellipse is P.A.Ri 63° (N over E). 



3.1. Substructure 

Figure H clearly indicates that the core of CL 0939-1-4713 
does not consist of a single spherically symmetric struc- 
ture, but is composed of two subclusters: Ml and M2. 
In Fig. 121 (left panel) a surface brightness projection of the 
cluster is shown, which is extracted from the region shown 
in the right panel. The two extended substructures form- 
ing the core are of similar size and of similar brightness, 
the subcluster in the East -M2- being slightly brighter 
than the subcluster in the West -Ml-). They both have 
an irregular shape, showing signs of further internal sub- 
structure. Their shapes are roughly elliptical with minor 
axis of « 120 kpc (19") and major axis of « 180 kpc 
(28").^ Their peaks are « 280 kpc apart and their coordi- 
nates (J2000) are: 

Pmi: RA= 9'^ 42"^ 58.9% DEC== +46° 59' 36.0" 



M2' 



RA= 9*^ 43°" 03.3% DEC= +46° 59' 26.5" 



They are slightly offsetted (Pmi ~ 44" north and 
Pm 9 ~ 21" south) comp ared to the position observed 
by ISchindler et al.l l)l998|) . and reach a maximum sur- 
face brightness of 5.0 x 10~^ counts s~^ arcmin"^ and 
5.8 X 10~^ counts s^-^ arcmin^^, respectively. 
The smaller peak in Fig. |21 (left panel), located on the 
left of the projection at a distance of about 90" from the 
left b order of the region, is caused by the presence of a 
QSO llDressler et al.lll993l) (see §E3|). 



3.2. Surface brightness profile 

Although the cluster morphology is not spherically sym- 
metric we perform a standard surface brightness analy- 
sis to obtain approximate values, which give a quantita- 
tive indication of the size of the cluster and also for com- 
parison with other clusters. We generate an azimuthally 
averaged surface brightness profile for the cluster in the 
0.3 — 2.0 keV energy band for each camera. This energy 
band was selected to optimize the signal— to— noise (S/N) 
ratio. Point sources were removed manually and the pho- 
tons were binned into concentric annuli with a width of 
3.3" (3 pixels of the MOS camera). The vignetting cor- 
rection and the background subtraction was done as de- 
scribed in § 12.3! and § 12.41 The three profiles were then 
added. 

As the cluster is clearly not spherically symmetric the pro- 
files were centred on the positions of the two peaks (Pmi 
and Pm2) of both subclusters. The profile centred on Pm2 
is obtained excluding a sector of 120° in the direction of 
the other subcluster. In the case of the profile centred 
on Pmi a wider sector was excluded (240°) in order to 
avoid the emission coming from the region located at the 
southern- west part of Ml. The resulting surface bright- 
ness profiles S{r) are shown in Fig. O where the overall 
XMM point spread function (PSF) binned as the observed 
profiles is plotted as a dashed line. The overall PSF is ob- 
tained by adding the PSF of each camera, estimated at 
an energy of 1 keV and weighted by the respective cluster 
count rate in the 0.3 — 2.0 keV energy band. The compar- 
ison of the surface brightness profile data with the PSF 
shows that the effect of the PSF on the profile is negligi- 
ble. 

After the subtraction of the corresponding blank field pro- 
file we find no residual CXB background left (as already 
mentioned in Sec. l2.4|l . We therefore fit the surface bright- 
ness with a /3-model without the need of an additional 
background component: 



S{r) = So 



1 



-1 -3/3+1/2 



1' = 389 kpc at the cluster distance of z = 0.406. 



where Sq is the central surface brightness, rc is the core 
radius, /3 is the slope parameter. The profiles for the two 
subclusters Ml and M2 are fit out to radii of 3.5' and 
3.2' respectively, at which the S/N > 3a. Results of the 
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Fig. 4. X-ray image of the galaxy cluster CL 0939+4713 in the energy range 0.3-2.0 keV. The image is a mosaic of the 
three MOSl, M0S2 & pn cameras. The core of CL 0939+4713 does not show a simple single central structure, but 
is composed of two subclusters. The two subclusters in turn show internal structure. Several other point sources are 
visible. The X-ray image is adaptively smoothed with a Gaussian with (Tmax — 6 arcsec to keep an uniform signal-to- 
noise ratio of 3 a. The contours are linearly spaced with a spacing of 0.006 counts s"^ arcmin^^. The highest contour 
line corresponding to 0.10 counts s~^ arcmin"^. 



different fits are given in Table ^ We obtain similar fit 
parameters for the two subclusters; however the central 
region of the profile centred on Pmi is not well fitted 
probably because of internal structure of the subcluster 
that flattens the central part. The derived beta and core 
radii parameters for both subclusters are larger compared 
with the ones obtained by Schindlcr ct al. (199§}. These 
discrepancies could be due to the low S/N ratio of the 
HRI observation that provides artificially low values for 
the beta-parameters. The HRI p rofiles decrease r a pidly to 
background level (see Fig. 2 in ISchindler et al.1 ljl998l) ). 
therefore a beta model fit is actually probing only the 



central part of the profile. In order to test this possi- 
bility decreasing cut off radii were applied to the XMM 
profiles. Smaller values of the core radius and beta were 
indeed derived: for a cut off radius of 2' Tc = 0.5' and 
(3 = 0.41 were obtained for the profile centred on Pmi and 
Tc — 0.8' and /3 = 0.57 for the profile centred on Pm2- 
For a more detailed analysis of the substructure, we sub- 
tract the emission coming from the subcluster M2 from 
the original image by using the /3-model with the fit pa- 
rameters in Tabled The residual image is shown in Fig. El 
In this residual image the emission from the quasar and 
from the point sources that appear in Fig.0|is still visible. 
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Fig. 5. Combined MOSl, M0S2 and pn surface brightness profile around the maxima M2 (left) and Ml (right) with 
spherically symmetric /3- model fits (solid line). The fit parameters are given in Tabled The dashed line shows the 
on-axis PSF, normalized to the central intensity. 

Table 1. Fit parameters for a (3 model centred on Ml and M2: 5*0 is the central surface brightness, Vc is the core 
radius, and /3 is the slope parameter. 



region 


in counts s ^ arcmin 


rc in arcmin 


13 


M2 
Ml 


8.6±0.4 X 10'^ 
6.8±0.5 X 10"^ 


1.11±0.07 
1 lfi+°-i* 

^■^"-0.15 


0.68±0.03 

n 66+°"^ 

^^•""-0.05 



There is clearly extended emission left from the subcluster 
Ml, but its peak has shifted to the south-west of about 
25 arcsec. North of Ml there is also residual emission sur- 
rounding some point sources, which is clearly extended 
and structured. South-west of Ml, separated by an arti- 
fact minimum caused by the gap in the detector chips is 
a last piece of residual extended emission. 

3.3. Count rate 

The X-ray emission of the cluster can be traced out to 
3.5 arcminutes (1.4 Mpc). The MOS and pn count rates 
within the R < 3.5' region and in the considered energy 
band (0.3-2.0 keV) are 0.125±0.013 counts s'^ and 0.31± 
0.03 counts s^^ respectively. 

3.4. Other sources in the pointing 

In the pointing many other sources are found around the 
cluster. 19 point sources, located in the area covered by 
MOS CCD#1, were detected by both the pipeline detec- 
tion and a peak density finding algorithm. Their coordi- 
nates, extension and count rates (after background sub- 
traction) are listed in Table |2| Most of these sources are 
shown in Fig. |3 (right panel) together with their identifi- 
cation number. 



For the sources corresponding to member galaxies, or lo- 
cated within 4.3 arcmin from the cluster centre, a more 
detailed analysis has been performed. 

ID #3 For this source no optical counterpart was found. 
The low statistics of the data do not allow us to con- 
strain the photon index, even using a simple power 
law model. Assuming a power law spectrum with pho- 
ton index F = 1.7 we obtain an unabsorbed flux of 
/x(0.3 - 10.0) = 5.7 X 10^15 gj.g s-icni-2. 

ID #4 No optical counterpart. Assuming a power law 
spectrum with photon index F — 1.7: /x(0.3 — 10.0) — 
3.5 X 10^^*^ erg s~^cm^^. 

ID 7^7 No optical counterpart. Assuming a power law 
spectrum with photon index F = 1.7: /x(0.3 — 10.0) — 
2.2 X 10"" erg s^^cm^^^ 

ID #8 Probable member galaxy Abell 0851: [DG92] 124. 
Assuming a power law spectrum with photon index 
F = 1.7: /x(0.3 - 10.0) = 1.0 x 10'^^ erg s-^cm-^, 
and a bolometric luminosity of _Lx(bol) ~ 1.6 x 
10"^ erg s-^ 

ID #9 Member galaxy Abell 0851:[DG92] 311 at a red- 
shift of z = 0.4007. It is also a radio source known as 
[CDS90] R082. Its spectrum was fitted with a power 
law giving a photon index F = 1.6 ± 0.2; the unab- 
sorbed flux is /x(0.3-10.0) = 4.7x10-1'^ erg s-^cm-^ 
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Table 2. List of point sources detected in the field covered by MOS CCD#1. Columns 1-3 give the source identification 
number, and its coordinates (J2000). In Cols. 4 and 5 the source count rate (in counts per second) and the radius (in 
arcseconds) out to which the source can be detected, respectively, are listed. 



ID 


RA 




DEC 






count rate 


radius 
["] 


1 


9" 


42"^ 


29.23" 


+46° 


56' 


25.27" 


0.0025 


30 


2 


9'' 


42'" 


34.90" 


+46° 


55' 


58.35" 


0.0024 


30 


3 


9'' 


42m 


40.53" 


+47° 


01' 


40.41" 


0.0006 


22 


4 


9*^ 


42'" 


41.72" 


+46° 


58' 


32.43" 


0.0004 


20 


5 


9" 


42'" 


44.04" 


+47° 


03' 


36.45" 


0.0016 


25 


6 


9'' 




51.09" 


+47° 


05' 


14.49" 


0.0028 


35 


7 


9" 


42'" 


51.88" 


+47° 


02' 


04.49" 


0.0023 


25 


8 


9'' 


42"' 


55.21" 


+46° 


57' 


10.50" 


0.0011 


17 


9 


9" 


42"' 


57.36" 


+46° 


58' 


46.50" 


0.0048 


15 


10 


9'' 


42"' 


57.39" 


+47° 


00' 


50.12" 


0.0054 


22 


11 


9" 


43m 


01.85" 


+47° 


01' 


24.50" 


0.0078 


30 


12 


9'^ 


43'" 


02.25" 


+47° 


02' 


16.50" 


0.0025 


22 


13 


9'' 


43'" 


03.03" 


+47° 


03' 


52.50" 


0.0026 


25 


14 


9'' 


43m 


04.20" 


+46° 


58' 


40.49" 


0.0036 


15 


15 


9'' 


43"' 


05.36" 


+46° 


56' 


36.49" 


0.0019 


18 


16 


9" 


43"' 


07.14" 


+47° 


02' 


16.48" 


0.0025 


25 


17 


9'' 


43" 


07.72" 


+46° 


59' 


40.48" 


0.0036 


15 


18 


9'' 


43m 


08.88" 


+46° 


57' 


00.45" 


0.0021 


27 


19 


9'' 


43"' 


17.69" 


+47° 


00' 


38.40" 


0.0089 


30 



and the bolonietric luminosity _Lx(bol) = 7.5 x 

ID #10 Probable member galaxy Abell 0851:[SED95] 
074. The spectral fit yields T = 1.3 ± 0.1; /x(0.3 - 
10.0) = 6.5 X 10"" erg s^^cm-^; ix(bol) = 1.1 x 
lO'^^ erg s-\ 

ID #11 This source was already detected during a 
ROSAT PSPC observation as RXJ0943.0+4701 
by .Schindler fc WambsganssI l)l996(l and had a mea- 
sured flux in the ROSAT (0.1 - 2.4 keV) band of 
/rosat,pspc = 4 X 10"^'' erg s'^cm-^. During 
a s econd ROSAT ob serva tion of the cluster done 
by ISchindler et"al] l(l998l) using the HRI, the 
source was hardly visible with a resulting flux in 
the (0.1 - 2.4 keV) band of /rosat.hri = 5 x 
10~^^ erg s~^ cm~^, indicating that this source has 
a very strong variability. The source coordinates ob- 
served by PSPC are shghtly offset (« 38" south) with 
respect to the coordinates observed with XMM. 
We find the galaxy Abell 0851:[SED95] 066 as a pos- 
sible optical counterpart, which is a probable member 
galaxy The spectral fit yields P = 1.8 ± 0.1; /x(0.3 - 
10.0) = 6.8 X 10-" erg s'^cm'^; Lx(bol) = 1.2 x 
10'*'' erg s"^ We compute a flux in the (0.1 - 2.4 keV) 
band of /xmm — 4.7 x 10^*^ erg s^^cm^^, which 
is only slightly higher than the values computed 
bv lSchindler fc Wambseansd l)l996|) in their first X-ray 
observation of the cluster. Considering the changing 
values of the luminosity and the resul ts of the spec- 
tral analysis, the conclusion drawn bv i Schindler et all 



l)l998j) . of source #11 being a strong X-ray variable 
AGN, is, up to now, the most reliable one. 
ID #12 No optical counterpart was found for this source. 
The spectral fit yields P = 1.95 ± 0.2; /x(0.3 - 10.0) = 

2.0 X 10"" erg s'^cm-^. 

ID #13 No optical counterpart. The spectral fit yields 
P = l-Ttof, /x(0.3-10.0) = 2.5x lO-i'^ erg s-^cm-^. 

ID #14 Probable member galaxy AbeU 0851: [DG92] 292. 
The spectral fit yields P = 1.6 ± 0.1; /x(0.3 - 10.0) = 
3.3 X 10"" erg s^^cm^^; Lx(bol) = 5.4 x lO*^ gj.g s"!. 

ID #15 Probable member galaxy AbeU 0851: [DG92] 064. 
Assuming a power law spectrum with photon index 
P = 1.7 we obtain an unabsorbed flux of /x(0.3 — 
10.0) — 1.8 X 10"*^ erg s~^cm~^, and a bolometric 
luminosity of Lx(bol) = 2.9 x 10^'^ erg s~^. 

ID #16 No optical counterpart was found. Assuming a 
power law spectrum with photon index P — 1.7: 
/x(0.3 - 10.0) = 3.7 X 10-" erg s-^cm-^. 

ID #17 This source is a QSO (known as AbeU 
0851: [DG92] 440) lying in the fie ld of view of the clus - 
ter and at a redshift of z = 2.055 ijPressler et al 
The spectral fit yields P = 1.72 + 0.05; /x(0.3-10.0) = 

3.1 X 10"" erg s'^cm'^; Lx(bol) = 2.3 x 10*^ erg s"!. 
In the ROSAT energy band the QSO has a luminos- 
ity of LxMM (0.1 - 2.4) 8.6 X lO"* e rg s-\ which is 
slightly lower than was measured by ISchindler et alJ 
t^l9m (Lrosat (0.1 - 2.4) = 1 4 X lO^'^ erg s" ^ ). Due 
to the low number of counts ISchindler et al.l l(l998^ 
could not fit any spectral model to the data and, to 
compute the luminosity, they fixed the photon index 
at P = 2.3. Using the same value for the photon index 
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Fig. 6. Residual image in the energy band 0.3-2.0 keV 
after subtracting the spherical /3 model centred on sub- 
cluster M2 from the original image. The positions of Ml 
(plus sign) and M2 (diamond point) are indicated. After 
the subtraction there is extended and structured emission 
around the centre Ml. The emission from the quasar and 
the point sources that appear in Fig.0]is visible. The con- 
tours are linearly spaced with a spacing of 0.005 counts 
s~^ arcmin"^. The highest contour line corresponding to 
0.10 counts s~ 



the two luminosities, computed by ROSAT and XMM, 
are perfectly consistent. 

ID #18 Probable member galaxy Abell 0851:[DG92] 105. 
Assuming a power law spectrum with photon index 
r = 1.7 we obtain an unabsorbed flux of /x(0.3 — 
10.0) = 1.9 X lO"^"* erg s~^cm~^, and a bolometric 
luminosity of Lx(bol) = 3.1 x 10^^ erg s^^. 

ID #19 J0943.2-h4700; has no optical counterpart. The 
spectral fit yields T = 1.8 ± 0.1; /x(0.3 - 10.0) = 7.8 x 
10"^'* erg s~^cm~^. 

3.5. Comparison with the optical waveband 

In Fig. [7| an HST image of the central area of the cluster 
is shown. On it, superposed, are the X-ray contours. 
An extensive redshift survey of the galaxie s in the field of 
view o f CL 0939-1-4713 was performed bv iDressler et alJ 
l)l999j) and at least 60 galaxies were spectroscopically 
proven to belong to the cluster. Among these galaxies, 
those falling within the field of view of HST, can be seen 
in Fig. d surrounded by rectangles. 

Some of these galaxies (e.g. #9) are very luminous also 
in X-rays but, nevertheless, the X-ray luminosity of the 
central region of the cluster is strongly dominated by 
extended hot emission from the intra-cluster gas, peaked 
in the two subclusters. 

Both the galaxy number density distribution and the 



X-ray emission show two main peaks. The peaks of 
these two distributions are in similar regions. It is also 
int eresting th a t the mass distribution (as computed 
bv ISeitz et al.l l|l996|) using a weak lensing analysis) has 
a similar bimodal structure. 

The subcluster Ml coincides with the strongest maxima 
in both the galaxy number density and the mass distri- 
bution. The peaks of the galaxy density and of the mass 
distributions are « 26 arcsec West and « 15 arcsec South 
of the X-ray peak Ml. 

The other subcluster, M2, is located near the second 
maxima in the galaxy number density and in the mass 
distribution (their peaks are, respectively, « 16 arcsec 
west and ~ 20 arcsec south-east of the X-ray one). 



4. Spectral analysis 

For the whole spectral analysis, the ready-made on-axis re- 
sponse files (ml_r5_alL15.rsp for MOSl, m2_r5_alL15.rsp 
for M0S2 ^ and epn_ff20^Y9.rmf for pn '') were used. ^ 
Source and background spectra were accumulated sep- 
arately for each detector. Spectra were re-binned in 
order to obtain S/N > 5 a in each bin after background 
subtraction. 

As the spectra are already corrected for vignetting, spec- 
tra of the same physical region, observed with the three 
different detectors, can simply be fitted simultaneously in 
order to maximize the signal to noise ratio. 
All the spectral analysis was performed in the 
0.3 — 10.0 keV band. Data below 0.3 keV have been 
excluded to avoid residual calibration problems in the 
MOS & pn response matrices below these energies. 



4.1. Results 

The spectra for the cluster analysis were extracted 
from a circular region centered at RA = 09*^ 43™ 01.0^ 
DEC = +46° 59' 37.5" (J2000) and of radius r=3.0'. The 
emission coming from all the point sources inside the 
region was excluded. 

X-ray temperatures, metal abundances, fluxes and 
luminosities, were computed with XSPEC fitting the 
d ata to an absorbed isot hermal plasma emission code 
Kaastra fc Mewj lll993l) including the FeL calculations 
Liedahl et all ^|l99,'^^ '^. folded through the propriate 



of 
of 

response matrices. In Table 01 is a detailed description of 
the results obtained fitting the theoretical model to data 
from the three detectors at the same time, and for each 
detector independently. 



Both MOS files are for patterns from to 12. 
pn response file is for single-pixel events only. 
^ All response files can be downloaded via anony- 
mous ftp at xmm.vilspa.esa.es in the directory 
/ pub/ccf/constituents / extras /responses. 
* MeKaL code from now on. 
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Fig. 7. HST image of the central region of CL 0939+4713 ijBuson et al.ll200(]|) . Superposed are the (0.3 - 10.0 keV) 
XMM X-ray contours linearly spaced from 0.24 to 0.98 of the maximum. White rectangles surround galaxies known 
to be cluster members. 
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Fig. 8. MOSl, M0S2 & pn spectra with MeKaL model 
superimposed (plotted in different shades of gray). 

Temperatures obtained using MOS and pn detectors 
independently are largely consistent. The fit results for 
the redshift a re consistent with it s spectroscopic value 
computed by IPressler et al.l l|l999() {z = 0.406); good 
agreement is also observed between the fit value of 
the hydrogen column density and its Galactic v alue of 
nn = 1.24 X 10^° cm'^ ijPickev fc Lockmanlll990|) . When 
fixing these two parameters the best fit temperature and 
metal abundance of the cluster based on a simultaneous 
fit of all (MOS+pn) the spectra (using independent 
normalizations for the MOS and the pn cameras) are 
kT = 4.5 ± 0.2 keV and Z = 0.20 ± 0.06 Zq, respectively 
(the corresponding MOSl, M0S2 & pn spectra with the 
fitted MeKaL model superimposed are plotted in Fig.lHJ. 
Using the best fit model parameters we derived, 
out to a radius of 210", an unabsorbed flux of 
/x(0.3 - 10.0) 1.3 ± 0.1 X 10-12 crgs-icm-2, 

corresponding to a luminosity, in the cluster rest frame, of 
Lx(0.3-10.0) = 1.1 ±0.1 X lO''^ erg s'^ and a bolometric 
luminosity of ix(bol) = 1.3 ± 0.1 x lO''^ erg s^^ 



4.2. Temperature map 

Since the cluster has been observed to have an irregular 
morphology, a deeper analysis of its internal structure 
and dynamical state is required in order to gain a better 
understanding of it. 

The cluster was subdivided into 18 regions; their sizes 
have been chosen such that each region (after background 
subtraction) has similar number of photon counts. We 
fixed the HI column density and the redshift to the 
average values for the whole cluster and then fitted the 
spectrum for each region to obtain the values of the 
temperature. 

The resulting temperature map is shown in Fig. El where 
different shades of gray refiect the different temperatures 
of the various regions. Superposed are linearly spaced 
X-ray contours. 

A clear trend of higher temperatures is observed in the 




Fig. 9. Cluster temperature map; superposed are hnearly 
spaced X-ray contours. Darker shades of gray indicate 
higher temperature. A trend of higher temperatures is ob- 
served in the central region, located between the two peaks 
(Pmi and Pm2) in the X-ray emission. The high tempera- 
ture region also extends in the NW and SE directions. 

region situated in the centre between the two clumps 
forming the core. The region of higher temperature also 
extends in the NW and SE directions. 
Average temperatures were also computed from spec- 
tra extracted from the central hotter region and 
from the two regions surrounding it. All three re- 
gions have elliptical shapes and the coordinates of 
their centres, and the major and minor axes are: 
M2. RA= 9'' 43°" 08.11^ DEC= +47° 00' 03.48" 
(J2000), a = 1.5', b = 1.2', P.A.= 140°; Central. 
RA= 9'' 43"" 01.58^ DEC:= +46° 59' 18.02" (J2000), 
a = 2.1', b = 0.4', P.A.= 159°; Ml. RA= 9'' 42"" 55.50", 
DEC= +46° 58' 35.50" (J2000), a = 1.4', 6 = 1.0', 
P.A.= 179°. A MeKaL model was used; HI column 
density and the redshift are fixed to the average values 
for the whole cluster. The resulting temperatures, which 
clearly confirm the trend observed in the temperature 
map, are kT = 'i-Stoi keV, kT = G.Otol keV and 
kT = 4:.7toi keV (for M2, central and Ml regions, 
respectively). Detailed results of the fits are presented in 
Table H 

The temperature map in Fig. |21 is similar to the predic- 
tions of simulations of cluster merging, for an almost 
head-on encounter of two subclusters of similar mass, 
short l y before the merging process fa icker & Sarazin, 
l200lt iRitchie & Thomasll2002t iTakizawal Il999l^ . Similar 
temperature structure has also bee n observed in clu sters 
such as A399 and A401 (Marke vitch et alJ \l99^ and 
A3528, A1750, and A3395 (iDonnellv et a].ll200lh . 
The central hotter region therefore delineates an area 
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Table 3. Results of the spectral fits for the whole cluster. Column 1: line number. Column 2: detectors in use. In Cols. 
3-6 the fit parameters (temperature in keV, metallicity in solar units, redshift and HI column density in 10^" cm~^) 
are listed. In Col. 7 is the fit reduced . 

Detector Temperature Metallicity Redshift HI column density x'^/d.o.f . 
kT in keV in solar units [lO^° cm~^] 



1 MOSl,MOS2 & pn 4.5 ± 0.2 0.20 ± 0.06 0.406 (fixed) 1.24 (fixed) 1.09 

2 MOSl,MOS2 & pn 4.6 ± 0.2 0.22 ± 0.06 0.418lo E;j^ 1.24 (fixed) 1.10 

3 MOSl,MOS2 & pn 0.24+° 0.411°-^'^ 0.3l°-^ 1.10 

4 MOSl,MOS2 4.6 ± 0.3 0.23 ± 0.08 0.406 (fixed) 1.24 (fixed) 1.11 

5 MOSl,MOS2 4.6 ±0.3 0.24 ± 0.08 0.41 ± 0.01 1.24 (fixed) 1.11 

6 pn 4.4 ±0.4 0.16t;5 [;^ 0.406 (fixed) 1.24 (fixed) 1.09 

7 pn 4.5 ±0.4 0.19 ± 0.1 0.43 ± 0.02 1.24 (fixed) 1.10 

8 pn 5.lto;3 0.21 ±0.1 0.43 ± 0.02 ^-^tW 1.10 



Table 4. Results of the spectral fits for the three subregions using a MeKaL model. HI column density and redshift are 
fixed to the average values for the whole cluster. Column (1): region analyzed. In Cols. (2) and (3) the fit parameters 
(temperature in keV and metal abundance, respectively) are listed. In Col. (4) is the reduced of the fit. 



Region 


Temperature 
kT in keV 


Metallicity 
in solar units 


xVd.o.f. 


M2 

Central 
Ml 




0.15 ±0.11 

'-'•■''5-0.12 


0.96 
1.03 
1.07 



where the two subclusters are compressing and heating 
intra-cluster gas, while moving one towards the other 
in the beginning of a merging process. Simulations 
by Schindler & Miiller (1993) show that this hotter 
region between the two merging structures begins to 
appear « 0.5 — 1.0 Gyr before the merging. Comparing 
this r esult with simulations performed bv Takizawa (I999|, 
l20onh allows a determination of the dynamical state: the 
two subclusters in CL 0939±4713 will collide in a few 
hundreds of Myr. 



4.3. Metallicity variation 

A higher metal abundance is observed in region Ml 
compared to the abundance in M2 and to the abundance 
in the central region (see Table 0}. While in Ml a 
metallicity of 0.33+q found, the metallicity in M2 
amounts only to 0.15 ± 0.11 in solar units. The errors are 
relatively large because only photons from small areas 
were used. The 90% errors overlap slightly. But there is 
a clear trend that the metallicity is not constant. The 
high-metallicity region Ml coincides with a strong peak 
in the galaxy number density distribution. It seems that 
the galaxies present in this region have enriched the 
intra-cluster gas with heavy elements. It is interesting 
that the high metal content is still localised within this 
region despite the active dynamical state of the cluster. 
In order to maintain such a high metallicity during the 
merger process, the subcluster Ml, on its way towards 
the subcluster M2, did obviously, not accumulate a lot of 



gas from the front and did not leave large quantities of its 
gas behind, always retaining the same gas. An additional 
explanation could be that the gas in this region has been 
metal enriched only recently, so that the motion of the 
gas has not had time yet to distribute the metals. 



5. Mass analysis 

The major difficulty to estimate the total mass in this 
cluster is its active dynamical estate. CL 0939±4713 is 
undergoing a major merger and the hypothesis of hydro- 
static equilibrium and spherical symmetry are in principle 
not valid. We therefore instead assume a bimodal model 
where each subcluster follows an isothermal beta-model. 
To estimate the gas and total mass of the cluster we anal- 
yse each subcluster separately assuming for each of them 
hydrostatic equilibrium and spherical symmetry. 
The integrated mass contained within a radius r for each 
subcluster can be calculated from the equation 

Mi<r)^-J^T,^J'-^^^^ + '^^ (1) 
/impG \ dmr amr / 

where pgas and Tgas are the density and the temperature 
of the intracluster gas, and r,fc,/i,TOp, and G are the ra- 
dius, the Boltzmann constant, the molecular weight, the 
proton mass, and the gravitational constant, respectively. 
To estimate the gas and total mass for each subcluster, 
the P parameters obtained in the fits of the profiles cen- 
tred on the peaks Pmi and Pm2 are used (see Table^l. We 
assume a constant temperature of 4.7 and 4.8 keV for the 
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subclusters Ml and M2, respectively. The two subclusters 
are very close each other: the outer radius of « 22" (140 
kpc) corresponds to half of the distance between the sub- 
clusters. Therefore we can compute the gas mass and total 
cluster mass out to this radius adding the masses corre- 
sponding to each subclustcr. Fig. llOlshows the integrated 
gas and total masses out to 140 kpc for the two subclus- 
ters independently and the numerical values are given in 
Table El 

For an estimation of the total cluster mass out to a 



Table 5. Numerical values for the mass determination in 
the two subclusters Ml and M2. 





gas mass (10^'^ Af©) 
within 140 kpc 


total mass M{r) (W'-'^Mq) 
within 140 kpc 


Ml 

M2 


0.56 
0.65 


4 2+U-« 

^•^-0.7 
A Q+O-B 




10'° 



0.06 0.08 0.1 
Radius (Mpc) 



larger radius we use the fit centred on the subclustcr M2, 
where the X-ray emission is higher. In Table the gas 
and total masses found at two different radii are listed. 
At radius rsoo = 0.85 Mpc the total cluster mass and gas 
mass are: Mtot (rsoo) = 2.3 x lO^^'M© and Mgas (rsoo) = 
0.49 X IO^'^Mq. The gas mass fraction tends to increase 
slightly outwards («17% at 0.5 Mpc and «22% at 1 Mpc). 
The tabulated errors on the total mass come from the un- 
certainty in the temperature and the errors in the /? pa- 
rameters. The emissivity of the gas in the XMM energy 
band 0.3 — 2.0 keV is almost constant inside the uncertain- 
ties for the estimate temperature. Therefore we can derive 
the gas density distribution without the uncertainty of the 
temperature estimate. 

The mass distribution of part of CL 0939+4713 has been 



Table 6. Numerical values for the mass determination in 
cluster CL 0939+4713 for two specific radii; the masses are 
computed using the fit parameters for the profile centred 
on subclustcr M2. 





mass within 


mass within 




0.5 Mpc 


1.0 Mpc 


gas mass (10^* M©) 


0.18 


0.67 


total mass M(r) (lO^M©) 


1.01 + 0.3 


2.98 + 0.7 



recons tructed with a weak lensing analysis bv lSeitz et al.l 
l)l996|) and a comparison with our cluster mass estimation 
would be interesting. Unfortunately the restricted region 
within which the total mass has been calculated with the 
gravitational lensing analyses makes a proper comparison 
impossible. 



Fig. 10. Integrated total mass (thick lines, top) and gas 
mass (thin lines, bottom) profiles for the subclusters Ml 
(solid fines) and M2 (dotted lines). 

6. Discussion of the merger hypothesis 

Several observed characteristics of CL 0939+4713, either 
in the X-rays and in the optical, point to a scenario in 
which the two subclusters Ml and M2, forming the clus- 
ter core, with masses of 4.2xlO^^M0 and 4.8xlO^^M0 re- 
spectively, are in the course of a merging process, moving 
one towards the other along a direction coincident with 
the cluster major axis. The cores of the subclusters will 
pass through each other in a few Myr. 

1) The spectral analysis of CL 0939+4713 has re- 
vealed an elongated region of heated gas located be- 
tween the two X-ray peaks, and almost perpendic- 
ular to the cluster major axis. Hydrodynamic sim- 
ulations ^chindler fc Mi illcr 1 9 §,^; . I^icker fc S arazia 
l200l[ iRitchie fc Thomajl2002t (Takizawalll999[) show 
that the gas located between two subclusters is com- 
pressed and heated during the approach of the subclus- 
ters. This hotter gas tends to expand in the direction 
perpendicular to the collision axis, therefore it shows 
an elongated shape perpendicular to the collision axis. 
The gas temperature and the elongation of the central 
region increase as the two subclusters get closer, while 
the bulk of the X-ray emission still comes from the 
cores of the two subclusters until the two cores collide. 

2) CL 0939+4713 shows an elongated X-ray morphology. 
Also this feat ure is consistent with simulations of merg- 
ing clusters llSchindleill2002j) . which show that before 
the merger the intracluster gas is elongated along the 
collision axis. 

3) The X-ray maxima of the two subclusters are displaced 
with respect to the peaks in the galaxy density distri- 
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bution. In Ml the X-ray peak seems to be ahead of 
the galaxy concentration. Obviously, the galaxies do 
not trace the dark matter distribution well in this re- 
gion. This is not surprising as during a merger the 
distributions of the different components can be quite 
different. 

4) CL 09394-4713 has a high velocity d ispersion, a k. 
1300 km s"^ l)Dressler fc GunnI Il992t jPressler et all 

indicating a relative motion of two or more 
components along the line of sight. The velocity dis- 
persions of the two subclusters, computed indepen- 
dently, are instead slightly lower: a ~ 1000 km s^^ 
and a ~ 750 km s^"'^ for regions M2 and Ml, respec- 
tively^; 

5) IPressler et al.l l|l999l) observed a very high fraction of 
post-starburst galaxies, in CL 0939-1-4713, compared 
to oth er rich clusters and to field galaxies. iFuiita et al., 
l)l999j) have investigated the effect of ram pressure in 
merging clusters. They f ind (in contrast to what was 
predicted previously by lEvrard 1,1991) ) that the in- 
crease in the star-formation rate caused by compres- 
sion is less important than the ram-pressure stripping, 
which dramatically increases before a merging process. 
The stripping of the interstellar medium causes a de- 
crease of the star formation rate, with a subsequent 
increase of the post s tarburst galaxy fraction. 

6) lAndreon et alJ ^^23) found that early-type galaxies in 

CL 0939-1-4713 do not show a radial symmetry, but are 
distributed along a privileged direction. This direction 
roughly coincides with the X-ray major axis and there- 
fore with the possible collision axis. 

7. Summary and Conclusions 

We confirm earlier results that CL 0939-1-4713 is a dv- 



namicallv voun} 


^ svste 


ISchindler et alJ 





The X-ray morphology has pro- 
nounced substructure. New findings of a hot region be- 
tween the two main clumps are particularly indicative of 
a merger process, in which the two major subclusters will 
collide in a few hundreds of Myr. 

Particularly interesting are the metallicity variations 
found in this cluster. Metallicity variations were found 



ies 1 


De Grandi & Molendil 


2001: 


Molendi fc GastaldcUd 


2001 


iKaastra et all l200ll 


Ga,sta,ldello fc Molendi. .20021 


Lewis et all l2002l). Chisters without cD galaxies 


were 



l200lt lArnaud et ani2001bj) . The problem when binning 
the cluster into radial bins is, that one averages over metal- 
licity values coming from very different regions of the clus- 
ters. Therefore it is not surprising that the profiles are 
usually flat for non-cD clusters. We took a different ap- 
proach here and found a hint that the metal abundances 
vary over the cluster. Of course it is quite difficult to do 
this in a relatively distant cluster. We therefore suggest to 
use this method in nearby clusters, instead of determining 
radial metallicity profiles. 
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